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The contribution of quarks to the effective potential for the phase of the Wilson 
line is computed at nonzero temperature and quark chemical potential to one and 
two loop order. At zero temperature, regardless of the value of the quark chemical 
potential, the effective potential for the phase of the Wilson line vanishes. At nonzero 
temperature, for special values of the phase the free energy of quarks equals that of 
bosons; at nonzero chemical potential, such quarks can "Bose condense", albeit with 
' negative density. 
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At nonzero temperature, SU(N) gauge theories in the absence of dynamical fermions 
have a deconfining phase transition related to the spontaneous breaking of a global Z(N) 
symmetry [|lj. The order parameter for this transition is the Wilson line: its vacuum ex- 
pectation value vanishes in the confined, low temperature phase, and is nonzero in the 
deconfined, high temperature phase. Because of the global Z(N) symmetry, in the decon- 
fined phase there are A" equivalent vacua, with the phase of the Wilson line equal to one of 
the N th roots of unity The interface tension between these A" vacua is computable 

semiclassically 0-[0], by computing the effective potential for the phase of the Wilson line. 



This interface tension is measurable from numerical simulations on the lattice flop . In the 
presence of dynamical quarks in the fundamental representation, the Wilson line acquires a 
vacuum expectation value at all temperatures, and the vacuum is unique, with zero phase 
for the Wilson line. Nevertheless, depending upon the number of colors and quark flavors, 
metastable states in which the phase of the Wilson line is nonzero can arise ||, and may be 
of physical interest (although this is controversial ||). 

The effective potential for the phase of the Wilson line for an SU(N) gauge theory has 
been computed at one and two PJ-P] loop order at a nonzero temperature, T. Here we 
extend this result to a nonzero quark chemical potential, \x. 
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The Wilson line is 



Q(x) = exp ig / A (x, t)cLt 



(1) 



where T is the temperature, g the gauge coupling constant, and the imaginary time r : — ► 
(3 = 1/T. To represent a nonzero phase of the Wilson line in the deconfmed phase, we can 
take a constant background value for the time-like component of the gauge potential, A 0- 
@ . After a global color rotation, this can be chosen to be a diagonal matrix: 
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The sum of the qj must vanish so that An, is an SU (N), and not a U(N), matrix. For example, 
consider Z(N) degenerate vacua || in the pure glue theory. In the trivial vacuum the phase 
of the Wilson line vanishes, f2 = 1, so all qj = 0. The equivalent vacuum in which the 
Wilson line is the first N th root of unity, fl = exp(2iri/N)l, is given by q± = . . . qjsr-i = q/N, 
and qN = — (N — l)q/N; tunneling between these two vacuum is parametrized by taking 
q : 0^ 1. 

Because a constant An. field can be eliminated by a global gauge transformation, the 
classical action is independent of the value of A . At nonzero temperature, this degener- 
acy is lifted by quantum effects, as the required gauge transformation alters the boundary 
conditions required for the fields in imaginary time. Consequently, at one loop order and 
beyond, a potential for Aq, or more properly for the phase of the Wilson line, is generated. 
In terms of partition functions ||, this potential is generated because one is summing over 
states which differ by a global gauge transformation. 

We compute the quark contribution to the effective action for a single, massless flavor. 
At one loop order, this is just the free energy in the presence of a background A field: 



S 1 = -Trlogi \l}> cl -//7o 



(3) 



Ij) cl = ql — igA 1 is the covariant derivative, and the trace includes color, Dirac indices, and 
a loop integral. We include the effects of a Fermi sea for quarks by introducing a chemical 
potential for baryon number, /i. 

In the following it helps to introduce two types of variables related to \i and the qj. One 
set of variables are color-dependent "chemical potentials" : 



Hj = n + litTi q^ , j = 1 . . . N . 

Notice that the qj act like an imaginary chemical potential for color 
complex, the effective action is as well. 

Alternately, we also use the variables Cj, 



(4) 

Since the \ij are 



C; 



1 . fi 

qj H h i^— 



J 



1...N . 



(5) 
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The fij are useful in comparing to the limit in which \i ^ and Aq — > 0; the Cj are convenient 
in the limit that A ^ and ji — > 0. 

As the background gauge field is diagonal in color, at one loop order the color trace is 
trivial, 

N N 

= 12 S J = - 12 lo S - Wo] • (6) 

We compute the Sj in two different ways. The first is careful and certainly correct, but 
awkward. The trace in Sj involves an integral over the loop momentum, = (k°,k), 
where for fermions in imaginary time, the energies are discrete, ko = (21 + 1)itT for integral 
I — ±1, ±3 .... Performing the sum over / by contour integration, 

Si = -PV 2T 3 J {log [l + e-( fc -^ T ] + log [l + e-^' T ] } ; (7) 

V is the volume of space, so (3V is the volume of space-time. To do the integral over spatial 
momentum, we use zeta function regularization; this was used at T ^ by Actor |TT|JT^| 
and by Weldon | 13[| . Expanding the logarithm and integrating over k gives 

Sj = -pv ^ £ ( -i)-«ijco S h (m) . ( 8) 

In each cosh function we expand as a power series, interchange the order of the two series, 
and so obtain a sum over rj functions, 

4T 4 » 1 /ft A2« 



where the rj function is defined as 

OO 1 

V(s) = j2(-l) i+1 - = (l-2 1 - s )as) Re(s)>l. (10) 

i=i 1 

Zeta function regularization is used to define the sum over negative s. Then a wonderful 
simplication enters: r](—2m) = for any integral m, so only three terms contribute to Sj, 

Using the explicit values of the rj function, 

( 2 2 TV j TV ) 

S^- P v\^N^ + T -^ i + ^-^Y (!2) 

For each this expression is periodic in qj modulo 1 over the interval G (—1/2, 1/2). The 
real part of S 1 agrees with Actor ||11||; he neglected the imaginary part. 
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The second way to compute S 1 is just to take the result for /i = 0, and substitute the 
Gj from d): 



4 I N 



(13) 



where BAx) is the Bernoulli function, 



B 4 (x) = x 2 (l - x) 2 ; 



(14) 



in this expression and henceforth, the real part of x is defined modulo 1. It can be shown 
that (|12|) and (p~3| ) agree. When the qj = 0, it is evident that (0) reduces to the free energy 
for a massless fermion at fi and [14j. When = 0, (|T3|) obviously reduces to the usual 
effective potential for the at T 7^ 0. 

While the first derivation was cumbersome, it has the virtue that in @ we only need the 
rj(s) function for negative values of s. In contrast, with the second derivation it is necessary 
to analytically continue the Bernoulli function BAx) for complex values of x, which is valid 
for \x — 1/2 1 < 1. Both (|l|) and (^) are analytic continuations of the same function, (f7|), 
at nonzero values of the chemical potential. 

Given our success at one loop order, however, we can follow the same process at two loop 
order. In other words, in order to compute the effective potential at two loop order, we just 
take the result for fi = and substitute the values of Cj in (||). This gives 



S 2 = - pv 9 —T^ 



N 



B 2 (C jk ) {B 2 {C 3 ) + B 2 (C k )) - B 2 (C J )B 2 (C k ) 



M _ 1 N /I 

■^ r T { {^-B2(C j ))B 2 (C j ) 



where we introduce the variable 



4 N 5 
tt E WO,) - B 3 (C k )) B x {C jk ) + ^(N 2 - 1) 

6 hk=l 



144 



(15) 



Cjk — Cj — C k 



(16) 



and the Bernoulli functions 



Bi(x) =x—~, B 2 (x) 



x -x + ~ , B 3 (x) 
6 



x 



2 X + 2 X] 



(17) 



it is useful to note that for real x, which is all that is needed here, B\( 
0, this reduces to the ~ g 2 free energy for /i,T^0 [|M| . 



-x) 



-BAx). When 



1j : 

At two loop order, quarks contribute to the effective potential from a single diagram in 
which a gluon crosses a quark loop. The first term in (|ToD, involving B 2 (Cjk), represents the 
case in which the two quark lines in the loop carry different Cj] then the exchanged gluon 
is off-diagonal, and Cj k = Cj — C k enters. (Notice that the quark chemical potential drops 
out of Cjk-) In contrast, the terms in ( ]T5| ) involving B 2 (Cj) arise when both quark lines 
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in the loop have the same Cj, so that the exchanged gluon is neutral with respect to the 
background field. This is true for the gluon in Feynman gauge. 

The origin of the terms ~ B\ and B 3 in flI5|) is due to two sources. The first is the 
longitudinal, and so gauge dependent, part of the gluon propagator, which contributes when 
the phase of the Wilson line is nonzero. 

The second source was noted by Belyaev || : what is physically relevant is not the effective 
potential for constant A$, but that for the phase of the Wilson line. This is most efficiently 
computed using a constrained functional integral approach, as discussed by Korthals Altes 
0. The effective potentials coincide at one loop order, but at two loop order, a finite 
renormalization of the Wilson line must be taken into account. This introduces again Bi and 
B3,: B3, for example, is the derivative of the one loop potential. The gauge dependence from 
this source cancels precisely against that of the longitudinal part of the gluon propagator. 

We draw two conclusions from our results. Consider first the limit of zero temperature. 
The parameters for the phase of the Wilson line, the qj, enters into the /ij ~ qjT, From 
the effective potential in terms of the /i/s, (|12"1), at one loop order the effective potential 
only depends upon the qj through terms which are at least quadratic in the temperature, 
~ ~ g|/x 2 T 2 . This remains true at two loop order, although it is less obvious for a 
potential written in terms of the Cj's, flT5|), than in terms of the fi/s. 

That the effective potential for the phase of the Wilson line is flat at zero temperature is 
reasonable. A phase for the Wilson line is parametrized by a constant A field; at nonzero 
temperature this cannot be undone by a gauge transformation, because it alters the boundary 
conditions at T 7^ 0. At zero temperature, however, the gauge transformation required to 
undo a constant Aq field is at infinite time, and so inconsequential. 

We stress that our result implies that only the phase of the Wilson line is independent 
of fi at T = 0. For a theory with dynamical quarks, the Wilson line has a nonzero vacuum 
expectation value at jj, — T — 0, with a phase which vanishes. When T = 0, as fi increases 
the vacuum expectation value of the Wilson line will as well, as the effects of baryon, or 
equivalently quark, loops enter. This increase is not calculable from the present work. What 
we do show is that nothing interesting happens with the phase. Remember that at T 7^ 
and /i = 0, there are metastable states associated with the phase of the Wilson line ||. Our 
results show that there is nothing analogous when T = and fi 7^ 0. 

Our second conclusion is when T ^ 0, any point at which a qj = ±1/2 is special. From 
the definition of (U), for such a qj the energy k° changes from an odd multiple of ttT (plus 
the chemical potential), into an even multiple. The corresponding part of the free energy, 
Sj, is then that of a free boson at |i,T/ 0. Thus whenever qj = ±1/2, the quarks "Bose 
condense": the true ground state is not a filled Fermi sea, as for any other value of qj, but 
a state in which there is a macroscopic density of states with zero momentum. Depending 
upon the q^s, several colors may condense at the same time. For two colors, for example, 
one can have all colors Bose condense at the same point. When qj = ±1/2, the two loop 
terms in the effective potential are just some perturbative corrections. 

This does not violate the spin-statistics theorem, because a system with qj 7^ does not 
satisfy standard thermodynamic properties. For example, from ((?]), the statistical distribu- 
tion function is n(k) = 1/ (exp((k — fij)/T) + 1); when qj = ±1/2, fij = // ± inT, and n{k) 
is minus the statistical distribution function for bosons! Because of behavior such as this, 
some authors have argued that qj 7^ is an artifact of the imaginary time approach |J . We 



5 



do not think so (see, e.g, Holland and Wiese, and Bronoff et al. 0), but do not presume to 
have ged the doubts of others. 

While we have only computed at one and two loop order, we believe that our two principal 
results — that the potential is flat for \i ^ and T = 0, and that there is "Bose condensation" 
for fij = ±1/2 — are true order by order in the loop expansion. In particular, in going from 
one to two loop order, there is no sign that the effective potential develops any singularity 
as T — > 0. 

The work of R.D.P. is supported in part by DOE grant DE-AC02-98CH10886. 
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